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ABSTRACT 

The standard picture of CV secular evolution predicts a spike in the CV distribution 
near the observed short-period cutoff Pq ~ 78 min, which is not observed. We show 
that an intrinsic spread in minimum ('bounce') periods P\, resulting from a genuine 
difference in some parameter controlling the evolution can remove the spike without 
smearing the sharpness of the cutoff. The most probable second parameter is different 
admixtures of magnetic stellar wind braking (at up to 5 times the GR rate) in a small 
tail of systems, perhaps implying that the donor magnetic field strength at formation 
is a second parameter specifying CV evolution. We suggest that magnetic braking 
resumes below the gap with a wide range, being well below the GR rate in most CVs, 
but significantly above it in a small tail. 
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1 INTRODUCTION 

The orbital period P is the one parameter of a cataclysmic 
variable (CV) which observers can usually measure with con- 
fidence. Roche geometry implies a close relation 
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between this period and the mass and radius A/2 , R2 of the 
mass-losing star. The accretion luminosity of CVs suggests 
that the mass transfer timescale tM = — A/2/A/2 is consid- 
erably shorter than the age of the Galaxy, so that M2 , P 
change on this timescale. Hence the observed distribution of 
CV periods is by far the most significant indicator of CV 
evolution. As is well known, it is consistent with the idea 
that the evolution is driven by angular momentum losses 
from the binary orbit (see e.g. King, 1988 for a review). 

The observed CV histogram (Fig. 1) cuts off sharply at 
an orbital period of P — Po — 78 min. There are respectively 
and 12 systems in the period ranges Po ± 5 min. The idea 
that Po represents a global period minimum (P = for Po) 
for CVs has been widely accepted for the last two decades. 
It is clear that such a global minimum can exist (Paczyhski, 
1981; Paczyhski & Sienkiewicz, 1981; Rappaport, Joss & 
Webbink, 1982). As the mass M2 of an unevolved secondary 
star in a CV is reduced by mass transfer, the binary period P 
usually decreases also, with R2 ~ M2 ~ P. However for very 
small A/2 <0.1Mq, the secondary's Kelvin-Helmholtz time 

iKH exceeds the timescale tM = —M2/M2 for mass trans- 
fer driven by gravitational radiation. Instead of shrinking 
smoothly to the main-sequence radius appropriate to its re- 
duced mass, the star cannot now reduce its entropy quickly 



enough and contracts more slowly. From (Jit) we see that 

1/3 

once R2 decreases more slowly than M 2 the orbital period 
P must begin to increase, defining a minimum ('bounce') 
period p, for the system. It is important to realize that 
there is nothing extreme about conditions in the binary at 
this point; in particular the mass transfer rate remains al- 
most precisely the same (~4x 10 _11 Mq yr _1 ) as when the 
system was well above the period minimum, cf Kolb (2002), 
Fig. 3. Nor do conditions change drastically thereafter: since 
M2 ~ O.IMq at Po we have tM ~ 2.5 Gyr. Hence most CVs 
do not have time to evolve to significantly longer periods or 
lower mass transfer rates after reaching Pb : clearly to evolve 
to —A/2 = 1 X 1O _11 M0 yr _1 would require of order the 
age of the Galaxy. CVs reaching Po thus remain clustered 
there with mass transfer rates similar to those of systems at 
longer periods. 

Detailed calculations always predict a value Pb very 
close to, if slightly shorter than, the observed Po. The dis- 
crepancy Pb < Po is persistent, but may reflect uncertain 
or over-simple input physics (cf Kolb & Baraffe, 1999). In 
particular the difference between the true and spherically 
approximated radii may account for most of the disagree- 
ment. However there is a much more serious problem with 
this interpretation of the observed cutoff at Po. This con- 
cerns the discovery probability 
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Here a is some (presumably positive) power describing ob- 
servational selection effects. For example a = 3/2 for a bolo- 
metric flux-limited sample, while a = 1 is often assumed 
for systems detected via optical outbursts. Since P = at 
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Figure 1. CV period histogram and cumulative distribution 
(dashed) below the observed period gap (data from Ritter & Kolb, 
1998). The Hcavisidc function shown is the representation (W). 
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Figure 2. Simple construction of a square distribution by adding 
discrete L-shaped individual components. The parameters used 
for the case shown are: r = 10, a = 0.83, AP = 1.5, Po = 78. 



P — P\,, p(P) must clearly have a significant maximum there 
unless — M2 declines very sharply near this period. In other 
words, the observed CV period histogram should show a 
sharp rise near a global minimum Pb unless the mass trans- 
fer rate drops there. However, as mentioned above, all evolu- 
tionary calculations show that — M2 changes very little as Pb 
is approached. We conclude that if Po is a global minimum 
there should be a large 'spike' in the CV period histogram 
there (cf Kolb & Baraffe, 1999). Instead, the observed period 
histogram (Fig. 1) has a 'corner' at Po, i.e. a sharp cutoff of 
a fairly flat distribution, but no spike. 

The lack of a spike in the observed distribution has 
prompted numerous theoretical investigations (see Kolb, 
2002 for a review, and Barker & Kolb 2002). Most of these 
propose ways in which CVs might become difficult to dis- 
cover near Po. A basic problem for this type of argument is 
that, as we have seen, there is nothing at all unusual about 
the system parameters (mass transfer rate, separation etc) 
at this period. Further, attempts to use accretion disc prop- 
erties as a way of making systems hard to discover founder 
on the fact that the AM Herculis systems, which have no ac- 
cretion discs, have precisely the same observed short-period 
cutoff Po ~ 78 min, and no spike either. King & Schenker 
(2002) suggested that CV formation may take roughly the 
age of the Galaxy, so that the oldest systems have not yet 
quite reached the minimum period. The fine-tuning here 
is perhaps worrying, although not easily disproved. In this 
paper we suggest a solution to the spike problem invoking 
a a different new ingredient in the standard picture of CV 
evolution. 



2 BLUNTING THE SPIKE 

The new element we introduce into the standard picture is 
the idea of a distribution of minimum periods Pb. At first 
sight this may not seem a particularly radical step. However 
it amounts to allowing the intervention of a second control- 
ling parameter. Second, note that the difficult part here is 



not smearing out the spike, but simultaneously retaining the 
sharp cutoff at Po . 

To fix ideas we first formulate a simplified version of 
the spike problem which has the virtue of being exactly sol- 
uble, and then proceed to a more realistic approach. In the 
simplified problem we make the following approximations: 

1. We assume that the relative discovery probability for 
any individual CV depends only on P — Pb, and can be 
represented as 



p(P - P h ) = H(P - P h ) - aH(P -P h - AP). 



(3) 



Here the TPs are Heaviside functions, Pb is the bounce pe- 
riod for the individual CV, and a, AP are constants. As 
can be seen by comparison with the lowest panel of Fig- 
ure 3 of Kolb (2002) this is a fair representation if we take 
a ~ 0.83, AP ~ 1.5 min, with Pb = 67 min in the case 
shown there. 

2. The observed CV histogram is taken as 



N(P) = N H(P - Po) 



(4) 



where No is a normalization. This is a relatively crude repre- 
sentation of the observed histogram (Fig. 1), but does cap- 
ture in essence its only significant feature, the 'corner' at 

Armed with these assumptions, we determine the rel- 
ative distribution n(Pb) of CVs with minimum period at a 
given Pb . Clearly the observed histogram N(P) is the con- 
volution of n(Pb) with p(P — Pb), so using (til ku) we have 



N(P) = / dP b n(P b )p(P - P b ) = N H (P - Po). (5) 
Jo 

Taking Laplace transforms, the convolution theorem gives 

(6) 



n(s)p(s) = — e sP ° 



Now since 
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Figure 3. Predicted distribution, assuming (a) that F\, varies 
over a 15 min range, starting from 78 min, with n(P b ) propor- 
tional to exp(— fc (Pb — Pd)/AP, k = 1.8, and (b) Pb varies over 
a 5 min range, with k = 0.6. In both cases, the intrinsic width of 
the spike is 1.5 min. 

we get 

AT „-sPn °° 

H ( S ) = I _ ae -sAP = N ° a e ■ (8) 

r = 

Hence transforming back we have 

oo 

n{P b )=N Y,a r S(Pb~Po-rAP). (9) 

r = 

This result gives a decomposition of the assumed histogram 
which is completely obvious when plotted (Fig. 2) . However 
it does illustrate the main point: one can create a rather 
flat but cut-off distribution (cf ^) from the very spiky indi- 
vidual probabilities (^) by taking a distribution of bounce 
periods P b > Po falling off with P b (here as a <Pb/AP> ) for 
Pb > Po- Dropping the higher-order terms r > r max pro- 
duces a downward step at P = Po + r m „AP in an oth- 
erwise uniform distribution, the discontinuity having size 
a rmax . If a 15% downwards glitch is regarded at the limit of 
acceptability when compared with observation, we can es- 
timate the spread in Pb required to wash out the spike as 
r max AP = (log 0.15/ log 0.83) x 1.5 min = 15 min. Similarly 
if we define ri/% by a 1 " 1 / 2 = 0.5 then more than 50% of CVs 
in the period range Po < P < Po + r^AP are close to 
their local minima P b (i.e. are in their local 'spikes'). We 
find ri/2 = 3.72, so we expect half of the CVs in a 5.6 min 
range above Po to be near their local P b values. 

It is now easy to guess a continuous form of n(P b ) which 
gives a reasonable representation of the corner as n(P b ) = 
exp[-0.124(P b -Po)]. The result is shown in Fig. 3, with the 
fitting extended to a maximum Pb equal to (a) Po + 5 min, 
and (b) Po + 15 min. As can be seen, extending to the longer 
period produces a much more acceptable representation. 

We can thus state our main result: the CV distribution 
near the short-period cutoff Po is well reproduced if CVs 
have a declining distribution of individual minimum periods 
Pb above the floor value Pq. 



3 THE HIDDEN PARAMETER 

In standard CV evolution (e.g. King, 1988) the global prop- 
erties of a system at any epoch are specified by two parame- 
ters, e.g. the initial primary and secondary masses Mu, Mn, 
corresponding to a white dwarf and an unevolved low-mass 
star respectively. Given these two quantities we can in prin- 
ciple calculate all others such as period, mean mass transfer 
rate and discovery probability P,—M-2,p(P) at any subse- 
quent time. In practice evolutions starting from different 
initial secondary masses M-a converge very rapidly to an ef- 
fectively common track (Stehle et al., 1996), and even the 
dependence on white dwarf mass is very weak. Thus in this 
standard picture all global properties of a CV are essentially 
specified if we know the orbital period P. This extreme sim- 
plicity gives the standard picture its predictive power, but 
also makes it vulnerable to problems such as the missing 
spike at Po. This is clearly predicted if every CV follows 
the common evolutionary track. To blunt the spike in the 
manner considered above requires sensitivity to a further 
physical parameter. Paczynski & Sienkiewicz (1983) show 
that this cannot be the primary mass or the assumed mean 
surface opacities (in any case the latter do not in principle 
constitute an independent parameter). 

We therefore turn to the other effect altering Pb consid- 
ered by Paczynski & Sienkiewicz (1983), namely increasing 
systemic angular momentum losses J by a factor / above the 
gravitational radiation value Jgr- They found a ~ 10 min 
spread in Pb for a range 1 < / < 2. We checked this using 
the codes of Mazzitelli (1989) - as adapted by Kolb & Ritter 
(1992) - and Hameury (1991), and found that the required 
> 15 min spread results for 1 < / < 6. 

Note that in fixing this range of / we assume that evolu- 
tion with f — 1, i.e. driven purely by GR, does produce the 
minimum bounce period Pb = Po — 78 min, and thus that 
theoretical efforts in this direction will succeed. We could 
instead in principle use a value / m i n > 1 to bring P b and Po 
into agreement, but this would rob us of the natural explana- 
tion for such a global minimum /mm, namely that this is the 
limit in which all other angular momentum loss mechanisms 
are small compared with GR. 



4 DISCUSSION 

We have shown that the corner in the CV period distribu- 
tion is reasonably well reproduced provided that CV evolu- 
tion has a second controlling parameter spreading the min- 
imum periods of an exponential tail of CVs by > 15 min. 
The most likely candidate for this parameter is an increase 
of the orbital angular momentum loss rate above that pro- 
vided by gravitational radiation by a factor /, ranging from 
1 to at least 6. Since gravitational radiation must always 
be present, the requirement / > 1 simply implies another 
angular momentum loss process whose strength ranges from 
values << | Jgr | to >5|Jgr|- An obvious candidate here is 
varying admixtures of magnetic stellar wind braking. The 
physical origin of the extra degree of freedom specified by 
/ might then simply be the intrinsic strength of the stellar 
magnetic field in the region of the ISM from which the sec- 
ondary star formed, which may then influence the strength 
of any dynamo-amplified or generated field. It is an ob- 
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Figure 4. Forbidden zone (shaded area) for angular momentum 
losses J in the interrupted magnetic braking picture (Rappaport, 
Verbunt &; Joss, 1983, Spruit & Patter, 1983) illustrated schemat- 
ically for a CV becoming fully convective at P ~ 3 hr with mag- 
netic braking rate Jo 3> 10 Jgr- To form the observed period gap 
between ~ 3 and ~ 2 hours, J must drop rapidly from a value 
of Jo to < 3 Jgri and stay low for a Kelvin time. After this, any 
amount of magnetic braking may be possible, leading to signifi- 
cantly different bounce periods f\, as given on the righthand axis. 



served fact that otherwise similar stars can have very dif- 
ferent magnetic field strengths, and this presumably applies 
to the secondary stars of CVs also. Apparently in most cases 
this field gives only weak magnetic braking [ Jmb] << | Jgr|, 
but there is evidently a tail of systems with stronger braking. 
This would for example naturally result if the distribution of 
intrinsic magnetic fields were a gaussian, with the mean suf- 
ficiently high to give strong braking above the period gap, 
but too low for such braking below the gap in all but an 
exponential tail of CVs. 

This is a satisfying conclusion, as one of the objections to 
the usual mechanism for forming the CV period gap be- 
tween 3-2 hours has been that strong magnetic activity 
is observed in main-sequence stars with masses low enough 
to be secondaries in CVs below the gap. The answer to this 
objection is now that indeed such activity may well occur 
in systems below the gap, but in most cases the resulting 
braking is weaker than GR. However for / ~ 6 much of the 
tail has values of magnetic braking considerably stronger 
than the extrapolation of the usual magnetic braking laws 
adopted for CVs above the gap (Verbunt & Zwaan, 1981 
Mestel & Spruit, 1987). Evidently these laws do not describe 
this regime well, in agreement with the idea that the nature 
of the braking changes character once the secondaries be- 
come fully convective (Taam & Spruit, 1989). It is therefore 
also plausible that intrinsically stronger fields are required 
below the period gap in order to give significant braking 
rates, thus accounting for the exponential tail of systems 
with strong braking as suggested above. 

We note that gap formation requires only that mag- 



netic braking Jmb should be interrupted, not suppressed, at 
P ~ 3 hr (plausibly where the secondary first becomes fully 
convective). Figure 4 illustrates cases where the braking af- 
ter the detached phase almost reaches the value before it. 
Provided magnetic braking does not resume, at something 
like its former strength, before the oversized secondary has 
managed to shrink significantly towards its main-sequence 
radius (on a timescale ~ £kh ~ 10 8 yr) , a gap of the observed 
properties will result. If the association with the change to 
a fully convective structure is maintained, we would thus re- 
quire that the resulting change in field topology or strength 
takes at least a time ~ £kh- To achieve the required smear- 
ing in Pb the angular momentum losses of a fraction of CVs 
must change in a way similar to the curve shown in Fig. 4. 

This discussion suggests that there is no obvious rea- 
son for any upper limit on the factor /, given only that the 
distribution of / values drops as / increases above unity. Sys- 
tems with large / have high mass transfer rates, and bounce 
at quite long periods. This is perfectly consistent with ob- 
servation, provided that only a minority of CVs below the 
gap have such high values of |Jmb|- We thus suggest that 
magnetic braking has the following properties (cf Fig. 4) 

(i) in most CVs above the period gap it drives mass 
transfer rates -M 2 >5 x 10~ 9 M Q yr" 1 

(ii) when a CV secondary becomes fully convective, 
magnetic braking rapidly drops to very low values for at 
least the thermal timescale of the star 

(iii) magnetic braking is present in all CVs below the 
period gap, and its strength depends on an inherent prop- 
erty of the secondary, possibly the formation magnetic field. 
In most cases it is weaker than gravitational radiation, but 
there is an exponential tail of systems with much stronger 
braking. These values are much bigger than given by ex- 
trapolating the usual magnetic braking laws to short orbital 
periods, and may thus require the strongest fields. 
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